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ABSTRACT 

We present a statistical parallax study of nearly 2,000 M subdwarfs with photometry and spec- 
troscopy from the Sloan Digital Sky Survey. Statistical parallax analysis yields the mean absolute 
magnitudes, mean velocities and velocity ellipsoids for homogenous samples of stars. We selected 
homogeneous groups of subdwarfs based on their photometric colors and spectral appearance. We 
examined the color-magnitude relations of low-mass subdwarfs and quantified their dependence on 
the newly-refined metallicity parameter, £. We also developed a photometric metallicity parameter, 
<5( g _ r ), based on the g — r and r — z colors of low-mass stars and used it to select stars with similar 
metailicities. The kinematics of low-mass subdwarfs as a function of color and metallicity were also 
examined and compared to main sequence M dwarfs. We find that the SDSS subdwarfs share similar 
kinematics to the inner halo and thick disk. The color-magnitude relations derived in this analysis 
will be a powerful tool for identifying and characterizing low-mass metal-poor subdwarfs in future 
surveys such as GAIA and LSST, making them important and plentiful tracers of the stellar halo. 



1. INTRODUCTION 

Modern s urveys, such as the Sloan Digital Sky Sur- 
vey (SDSS; I York et alj|2000l) and the Two Micron All- 
Sky Survey f2MASS: ISkrutskie et all l200fih . have pro- 
vided a detailed picture of the Milky Way's struc- 
ture, kinematics and evolution. The Galaxy's disks 
and halo have been measured by a slew of stellar trac- 
ers, including RR Lyrae stars, blue horizo ntal branch 
stars and main-sequence turnoff stars (i.e ., I Vivas et al.l 
[20Q]1 iBell et al.l 12001 lYannv et al.l I2009I ). But these 
tracers only represent a small fraction of the stars in 
the Milky Way. The most dominant stellar inhabi- 
tant of the Milky W ay, in mass and number, are the 
low-mass M dwarfs dReid et al.1 120021 : iChabrierl 120031 : 
iBochanski et al.l 120101 ) . Recent surveys have observed 
these stars in large numbers with deep, precise multi- 
band photometry over large fractions of the sky. The 
resulting photometric and spectroscop ic databases are 
the largest ever assemble d for M dwarfs ()Bochanski et al.l 
120101 : iWest et"al1l2011bj ) and have led to novel investiga- 
tions of the Galaxy. Most notably, the M dwarf stellar 
densi ty distributions (jBochanski et al.l [20101: Uur c et all 
200 8j ) and kinematic structures (IBochanski et a" 1 12007a!: 



Bond et all 120101 : iPineda et al.ll2012D have been used as 



tracers of the Galaxy's gravitational potential. 

The oldest members of the low-mass stellar popula - 
tions are the subdwarfs. First coined by Kuiperl (|1939l ). 
subdwarfs are low-metallicity stars ([Fe/H] < — 1) from 
earlier generations of star formation. The decreased 
opacity from metals leads to smaller stellar radii at the 
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same mass, as first noted by Sandagc & Eggcn ( 1959). In 
observational terms, low-metallicity subdwarfs lie below 
the main sequence of stars, with a fainter absolute magni- 
tude at the same color. Detailed studies of the low-mass 
subdwarfs are hindered by their low luminosities, but 
observational efforts have identified low-mas s subdwarfs 
throughout the M and L spectral sequence (iGizisl 119971: 
Burgasser et~aTI 120031: iDigbv et all 120031: lLodieu et ail 
2012). Subdwarfs are relatively rare in the solar neigh- 
borhood, with only 3 known wi thin 10 pc, compare d 
to 243 main sequence systems (|Monteiro et all [2006) . 
Traditionally, they have been identified by their large 
proper motio n s, such as Kapetyn's s t ar (u = 8.66"/ yr" 1 ; 
Mould! [19761: IWoolf fc Wallersteinl 12001 Ivan Leeuwenl 
20071) . Despite their paucity in the solar neighborhood, 



they probably comprise the majority of stars in the 
Milky Way's stellar halo. Precise metailicities for most 
M subdwarfs have not bee n measured, due t o the lack 
of adequate stellar models (jWest et al.ll2011al ). Instead, 
considerable effort has gone into developin g proxies for 
metal licity, most nota bly C Introduced b y [ Lepine et all 
(|2007ft and revised bv lDhital et~aTI (|2012l ). C is based on 
the molecul ar bandhead indices originally employed by 
Gizisl (|1997[ ) and crudely traces metallicity (Wo olf et all 
20091 ). While C does well at estimating [Fe/H] at low 
metallicity, it has limited precision near solar values. 
New empirical calibrations in the infrare d and optical 
(jRoias-Avala et al.ll2012HMann et aLll2012D may improve 
the calibration between £ and an absolute metallicity 
scale ([Fe/H]). 

Despite the advances in survey technology, precise dis- 
tances of SDSS low-mass dwarfs remain elusive. Most 
color-absolute magnitude relations (CMRs) used to ob- 
tain photometric parallax distances, are populated by 
~ 100 nearby stars, biasing the metalli city distribu- 
tion to values similar to that of the Sun dMonet et all 
fl^lRekTet al.lll99alBochanskil[20l)ll.Tao et alJ l2011). 
This is particularly troubling for calibrating the SDSS 
CMRs, with low-mass dwarfs detected out to a few kpc 
from the Sun. Co mplementary investigations of solar- 
mass stars in SDSS (jlvezic et al.ll2008HBond et al.ll2010D 
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have shown a gradient in metallicity with position in the 
Galax y, as well as a sh ift in the color-absolute magnitude 
locus ([An et all 120081 ). Using p hotometry and trigono - 
metric parallaxes of nearby stars. IBochanski et al.1 ()2010l ) 
derived color-absolute magnitude relations for M dwarfs 
in the uqri z filter set. However, n early all of the stars 
used for the IBochanski et al.l (|2010l ) relations are within 
20 pc and have spectra consistent with solar-mctallicity 
M dwarfs, limiting their applicability to low-metallicity 
subdwarfs. 

Usin g a statistical parallax technique. FBochanski et al.l 
(|201 If) estimated the absolute magnitudes for SDSS 
M dwarfs in the Data Release 7 catalog (DR7; 
lAbazaiian et al.l 120091: IWest et al.ll2011b[ ). They demon- 
strated that nearby low-mass stars are brighter in M r 
than most early-type M dwarfs in SDSS. This was at- 
tributed to the influence of metallicity and magnetic ac- 
tivity on the brightness of SDSS M dwarfs. However, 
the metallicity ra nge probed i n the s ample was relatively 
limited, since the lWest et al.l ()2011bf ) sample mostly con- 
tained solar-metallicity M dwarfs. 

In this paper, we present a statistical parallax analysis 
of over 2000 low-mass subdwarfs identified in the SDSS 
DR7 catalog. A summary of the observations employed 
in our analysis is contained in Section [5J The statis- 
tical parallax method, which yields absolute magnitudes 
and kinematics for a homogeneous groups of stars, is dis- 
cussed in Section [3] We find that subdwarfs with similar 
metallicity, traced by either spectroscopic or photometric 
indicators, form distinct color-magnitude relations, well- 
separated from the solar metallicity CMRs previously de- 
rived. We quantify these relations and compare subdwarf 
kinematics to main-sequence M dwarf kinematics in Sec- 
tion HI followed by our conclusions in Section [5] 

2. OBSERVATIONS 

The statistical parallax method relies on accurate and 
precise astrometry, photometry, proper motions and ra- 
dial velocities. All of these quantities are available within 
SDSS catalogs. Multi-band photometry was collected 
using the 2.5m SPSS telesc ope at Apache Point Ob- 
servatory (|Gunn et al.l I2006T) . The telescope employed 
a drift-scan technique, simultaneous ly imaging the sk y 
in ugriz along 5 camera columns (Gunn etaLl fT998h . 
Typical faint limits for the 53.9 second exposures were 
r ~ 22 mag, with bright objects saturating near r ~ 15 
mag. Typical systemat ic photometric precision is 0.02 
mag (jlvezic et al.l 120071 ) and has been quantified using 
repeat scans of "Stripe 82", a 300 sq. deg. equatorial 
patch of sky imaged repeatedly over the survey lifetime. 

Proper motions were supplied from the ProperMo- 
tions table contained in th e Catalog Archive Server 
(CAS: lO'Mullane et ai]|2005D . Proper motions were de- 
rived by matching SDSS detections to POSS observa- 
tions, with a typic al temporal baseline of ~ 50 years 
(jMunn et al.l 120041 ). This led to a precision of - 2 
mas yr _1 . Several cuts on precision and error flags 
were employed to ensure high quality proper motions. 
These are de scribed in detail in our previous investi- 
gations (i.e.. IDhital et al.l 120101: IBochanski et al.l 120111: 
IWest et ail 1201 lbTT Astrometric positions are reported 
by the SDSS pipeline with an internal precision of 25 
mas and an ab solute accuracy o f 45 mas in each coordi- 
nate direction (|Pier et al.1 12003'). 



When the sky conditions at Apache Point were not 
photometric, the SDSS telescope obtained R ~ 1, 800 
spectra of objects chosen by targeting algorithms de- 
signed to select hig h priority objects. This includes high- 
redsh ift quasars ([Richards et al. l | 2009t ISchneider et all 
120101) , luminous red galaxies (jEisenstein et all 120011) 
and o ther exotic target s, such as cataclysmic vari- 
ables dSzkodv et all 201 ID . The fiber- fed spectrographs 
(jUomoto et al.l 119*99 ) acquired 640 spectra simultane- 
ously. Over 600,000 stellar spectra have b een identi- 
fied i n the latest SDSS Data Release (DR9; lAhn et all 
120121) . We measured Ha emission and mol ecular band- 
head strengths using the Hammer pipeline (|Covev et all 
I2007D on the SDSS subdwarf spectra. Radial velocities, 
derived from comparing to the SDSS M dwarf templates 
(Boc hanski et al.ll2 007b) were measured for all objects in 
the lWest et all (|2011bl ) catalog, with a precision of ~ 7 
km s _1 . 

The subdwarfs in this an alysis passed the init ial color- 
cuts used to construct the IWest et all (|2011bl ) catalog. 
Subdwarfs were manually excluded from that catalog 
and placed into a separate list for further examination. 
There were 6600 spectra that were classified as "Odd" 
during the original spectral typing. These spectra were 
re-examined, and 260 low-mass subdwarfs were identi- 
fied by eye. We supplemented th is list by selecting 1767 
objects in thelWest et all (1201 lbl) catal og with ( < 0.825 
(|Lepine et all 120031: IDhital et al.l l2012f ). Further details 
and analysis of this subdwarf catalog can be found in 
iSavcheva etHI (|2013l ). 

3. METHOD 

The maximum likelihood application of the statistical 
pa rallax techniqu es employed here was first described 
in 
in 



Mur ray! (1 19831) and wa s applied to SDSS M dwarfs 
Bochans ki et all (|2011l ). It has also been used for 



previous studies of RR Lyrae stars (IHawlev et all 119861: 
I Strugnell et~aH Il986t ILavden et all 119961: iFernlev et al 
' 1998! IKol lmeier et all 120121), Cepheids (IWilson et al 



119911) and nearb y M dwarfs (IHawlev et al.l ll996ll We re- 
fer the reader to IBochanski et al.l (|2011[) for detailed de- 
scription of the method and briefly describe it here. The 
statistical parallax method models the kinematic prop- 
erties of a homogeneous population of stars with nine 
parameters, determining the mean velocities, dispersions 
and the orientation of the velocity ellipsoid. Two addi- 
tional parameters model the value and dispersion of the 
absolute magnitude of the population. The eleven pa- 
rameters used in the model are given in Table [1] In par- 
ticular, the distance scale and absolute ma gnitude values 
are gi ven by the following equation from IHawlev et all 

Gisg): 

a M 2 =log 10 [l + a fe 2 /(l + fc) 2 ]/(0.04 In 10), (1) 

where &m is the spread in absolute magnitude for a given 
luminosity bin, and is the spread in k which is the dis- 
tance scale parameter. The absolute magnitude is related 
to k through the following equation: 



M = 5 log 10 (l + k) + M A - 0.1 In 10(ctm) 2 



(2) 



where Ma is an initial estimate of the absolute magni- 
tude. 
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TABLE 1 

Statistical Parallax Fit Parameters 



Parameter 


Units 


Description 


<W 


km s — 1 


Velocity dispersion in radial direction 


a v 


km s — 1 


Velocity dispersion in orbital direction 


cr w 


km s — 1 


Velocity dispersion in vertical direction 


r, (f>, z 


radians 


Orientation of velocity ellipsoid 


U 


km s~ 1 


Solar peculiar motion (r) 


V 


km s~ 1 


Solar peculiar motion (</>) 


w 


km s — 1 


Solar peculiar motion (2) 


k 




Distance scale 


o"fc 




Dispersion in distance scale 



The positions, proper motions, radial velocities and ap- 
parent magnitudes, after corr ecting for Galactic ex tinc- 
tion using the dust maps of ISchleeel etall (1998), are 
input for each sample. The median r-band extinction for 
the entire sample was 0.07 mag, with 96% of the sam- 
ple having A r < 0.4. While most of our stars (99%) are 
more than 15 degrees away from the Plane, we compared 
A r for stars cl o ser to the Plane to the reddening maps 
of Uones et al.l (|2011[) . The m edian difference b etween 
our adopted reddening and the Uones et al.l 1)201 ID maps 
for these stars was 0.06 mag, well within the systematic 
uncertainties of this analysis. Observational uncertain- 
ties are used to properly weight the solution. The eleven 
parameters in the model were solved simultaneously by 
maximizin g the likelihood with a geometric sim plex opti- 
mization (|Nelder fc Meadll965tlDanielslll978l ). while the 
uncertainties were estimated using computed derivatives. 
The maximum likelihood e quations and s i mplex method 
are described in detai l in lHawlev et al.l (|1986l ). As in 
iBochanski et al.l (|201lD . we fixed <7fc at four values: 0.05, 
0.1, 0.2 and 0.3. We adopt at — 0.2 as a fiducial value, 
which corresponds to <tm ~ 0.4 which is the typical dis- 
persion in absolute magnitude for low-mas s stars in the 
solar neighborhood (Bochan ski et al.ll2010l ). 



3.1. Constructing Subsamples for Analysis 

Statistical parallax analysis requires the selection of a 
homogenous set of stars. Each set must be moderately 
populated (i.e., N > 30) to ensure a well-measured ve- 
locity ellipsoid. We performed a series of cuts on our 
subdwarf sample to explore the photometric and kine- 
matic properties of SDSS low-mass subdwarfs. They are 
listed in Table [H and described here. First, we sepa- 
rated the subdwarfs by r — z color. This color h as been 
shown to track M r well (|Bochanski et alj|2010j) . How- 
ever, metallicity can also influence the color-absolute 
magnitude loci. Thus, we further divided the sam- 
ple by r — z color and subdwarf subclass. The sub- 
dwarfs (sdM), extreme subdwarfs (esdM) and ultra- 
subdwarfs (u s dM) were chosen using the definitions of 
iLepine et all (|2007l ) and iDhital et ali (|2012f ). Example 
spectra of the different metallicity subclasses and their 
correspondence to C are shown in Figure [TJ The classes 
are separated by £, which is defined by the CaH and 
TiO optical bandhead strengths. Past studies have em- 
ployed the Ti O bandheads as tracers of absolute mag- 
nitude (i.e., iReid et all I1995D . while combinations of 
CaH and TiO have been used to trace metallicity shifts 



(lGizidll997t iBurgasser fc Kirkpatrick 2006; L epine et all 
l2007tlBochans ki et al.ll2007al ). As seen in Figure [TJ both 
CaH and TiO display sensitivity to metallicity differ- 
ence, with CaH near 6200 A and 6800 A strengthen- 
ing with decreasing metallicity, while the TiO bandheads 
weaken. Si nce low-metallicity su bdwarfs trace an older 
population (jMonteiro et al.ll2006[ ) , there is a considerable 
dearth of chromospheric Ha emission. Only one color bin 
(r — z = 1.65) contained enough active stars for analy- 
sis. Analysis of this bin, as described below , follow ed the 
general trend observed by IBochanski et ail (|2011f ). with 
active stars being brighter at the same color. 

While dividing low-mass subdwarfs into metallic- 
ity classes is useful for spectroscopic investigations, 
many current and planned large surveys, such as the 
Panoramic Sur yey Telescope fc Ra pid Response System 
(PanSTARRS; Kaiser et ail I2002D and Lar ge Synoptic 
Survey Telescope CLS ST; llvezic et al.l[2008D will consist 
solely of photometry. lHartwickl (|1977 ) originally noted 
that lower metallicity led to increased hydride absorp- 
tion in the B band giy i ng a redder B — V color for 
subdwarfs. iWest et ail (|2004f) demonstrated that the 
SDSS a — r color also re ddens with decreasing metal- 
licity. IWest et ail (|2011b|) further quantified this effect 
by fitting a two-dimensional relation between £, g — r 
and r — z (their Equation 3), but their relations are lim- 
ited to solar-metallicity M dwarfs. To quantify the shift 
in g — r for subdwarfs, we introduce S^ g ^ r j, a quantity 
that measures the difference in g — r between a subdwarf 
and its solar-metallcity counterpart as a function of r — z 
color. It is defined as: 

^(g-r) = (.9 _ r ) subdwarf - (.9 _ r) d M, (3) 

where (g - r) dM is: 



(.9 - r) d M 



6 

E 

i=0 



C n x(r- z) n , 



(4) 



with the polynomial coefficients Cq...Cq given in Table [3] 
and (g — r) suo dwarf is the g—r color of the subdwarf. The 
polynomial fit described in Equation|4]is shown in Figure 
[5] and was fit to the g — r,r — z color-color locus of the 
main sequence M dwarfs from SDSS DR9. These main 
sequence stars were selected to have good griz photome- 
try (i.e., 15 < griz < 20, a griz < 0.05, and Clean = 1). 
Carbon stars were removed by enforcing q — r < 1.8 for 
the fit (jMargon et al.l 120021 : ICovev et al.ll2007l ). Median 
g — r and r — z colors were computed, and a sixth order 
polynomial was fit from 0.8<r — z<3.6. 



The S 



metric is a crude measure of metallicity, 



and as discussed below, traces separate loci in color- 
magnitude space. Using the statistical parallax cali- 
brations discussed below, we will show that with grz 
photometry alone, one can estimate the metal content 
and distance of a low-mass star. Similar relations 
in other filter sets h ave been derived for nearby stars 
(|Johnson et al.ll2012l ). 

Finally, we removed RV outliers that were found more 
than 3 standard deviations away from the median RV 
of each subsample. Typically, this removed only one or 
two stars with large RVs ( > 300 km s _1 ). These were 
removed to prevent them from significantly influencing 
the mean velocity measured by the statistical parallax 
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TABLE 2 




SDSS SUBSAMPLES 




Initial Cut Bin Size Further Criteria 


Comments 


r — z 0.3 mag bins 




c 


sdM, esdM, usdM 


5 (9-r) 


0.1 mag bins 



Note. — Bins with fewer than 30 stars were not used in this 
analysis. 



TABLE 3 

(g — r) dM Coefficients 



Coefficient 


Value 


Co 


-14.8952 ± 3.0184 


Ci 


47.8830 ± 10.1991 


C-2 


-56.3038 ± 13.6807 


c 3 


34.1232 ± 9.35810 




-11.2926 ± 3.45636 


Cs 


1.94407 ± 0.656087 


c 6 


-0.136479 ± 0.0501870 



Note. — Polynomial fit to the 
g — r,r — z color-color locus of SDSS 
M dwarfs. Applicable to SDSS stars 
with 0.8 < r - z < 3.6. 



routine, since each RV was assumed to have a fixed un- 
certainty. F urther details on high velocity subdwarfs are 
discussed in lSavcheva et al.l ([2013D . 



3.2. Computational Method 

Using the method described in Bochanski et al. (2011) 
the following analysis was applied to each subsample. 
Ten loops were computed for each dataset, with 5,000 
simplex optimization iterations per loop, for a total of 
50,000 iterations per subsample. The initial absolute 
magnitude estimates for each sample were com puted 
from the M r , r-z CMR of lBochanski et al.l (|2010D . After 
each loop, the absolute magnitude estimate was updated, 
and the output of the previous run was used as input for 
the next. Typically, convergence was obtained after 5 
loops (i.e., 25,000 iterations). 

4. RESULTS 

There are three major results from our statistical par- 
allax analysis of subdwarfs: color-absolute magnitude 
relations, velocity dispersions and mean motions with re- 
spect to the Sun for each subsa mple. We compare our 
results to previous studies (i.e., iBochanski et al.l 120111) 
and discuss their relevance to large surveys of low-mass 
stars. 

4.1. Absolute Magnitudes 

For stars lacking trigonometric parallax measurements, 
photometric colors or spectral types have traditionally 
been employed to estimate their absolute magnitudes, 
resulting in photometric or spectroscopic parallax rela- 
tions. These relations, s uch a s those in lHawlev et all 
([200^1 . Uuric et al.l (|2008l) and IBochanski et al.l (|2010r ) 



are often based on the colors or spectra of nearby stars 
with trigonometric distance estimates. Deep surveys, 
such as SDSS, detect low-mass stars at distances out 
to a few kpc, which may not share the same age or 
chemi cal abundances a s M dw arfs in the solar neighbor- 
hood. Bocha nski et al.l ([201 ID demonstrated that age (as 
traced by Ha chromospheric emission) and metallicity 
(as traced by C) can produce differences of ~ f mag in 
M r at a given r — z color or spectral type. 

In Figure [31 we compare the st atistical parallax results 
from the current analysis and ( Bochanski et al.l 1201 1L 
green line) to the ne arby star photom e try. T he M r , r — z 
CMR measured by IBochanski et al.l (|2010t ) for nearby, 
solar-metallicity stars is shown with the red line. In- 
cluded in Figure[3]is the subdwarf LP 21-194, which was 
obser ved as part of the nearby star sample ([Bochanskil 
2 008 ). The subdwarf sample from the current study 
(blue line) shows a systematic shift compared to the 
main-sequence M dwarf statistical p arallax results and 
the ne arby stars. As first noted by iSandage fc Eggenl 
(1959), low-metallicity subdwarfs lie below the main se- 
quence of stars, with a fainter absolute magnitude at the 
same color. The subdwa rf sample does not span a com- 
parable color range as the IBochanski et al~l ([20111 ) results, 
since there are fewer stars to populate each color bin. 

To further quantify the dependence of absolute mag- 
nitude on metallicity, we calculated statistical parallaxes 
in the r, i and z bands. Both £ and ^( 9 -r) were used to 
trace metallicity, and the results are discussed below. 

4.1.1. Absolute Magnitude Variations and £ 

The left panels of Figures |H [5j and [5] show the de- 
rived CMRs for each metallicity subclass (sdM, esdM 
and usdM). The three classes are compared to nearby 
M dwarfs and the b est-fit nearby star CMR from 
IBochanski et al.l (|2010l ). The CMR for each subclass is 
roughly parallel to the main sequence, with an offset of 
~ 1 mag in absolute magnitude between adjacent sub- 
classes at the same r — z color. The largest separation 
from the main sequence occurs for the most metal-poor 
subclass, usdM. For example, at r — z ~ 1.35, there 
are ~ 3 mags difference between the main sequence and 
the usdM CMR. At a fiducial apparent magnitude of 
r ~ f8, this can result in overestimating the distance 
by a factor of ~ 5x. Thus, identifying the metallic- 
ity subclass is crucial for determining precise stellar dis- 
tance estimates. Despite the large spectroscopic cam- 
paigns p lanned for the near future such as GAIA and 
4MOST (|Perrvman et al.ll200ll: Ide Jong et al.lf2012h . the 
vast majority of subdwarfs will only have ugriz photom- 
etry, courtesy of LSST. 

4.1.2. Absolute Magnitude Variations and <5( 9 _ r ) 

A photometric proxy of metallicity for M dwarfs is 
complicated by the lack of a comprehensive spectroscopic 
calibration. However, the outlook is promising, as re- 
cent investigations have resulted in spectroscopic metal- 
licities of low-mass stars with ^0.1 dex precision, over 
limited rang es. These studies have focused mainly on 
the infrared ([Roias-Avala et al1l2012D. bu t some optical 
metallicity calibrations ([Mann et al.ll2012l) have been de- 
veloped. The ^(g- r ) metric was designed to use the ugriz 
photometric system. We eschewed the use of it, as the 
SED of M dwarfs is peaked near the near-infrared, with 
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dM C> 0.825 

sdM 0.5 <^< 0.825 
esdM 0.2<^<0.5 
usdM ^< 0.2 



4000 5000 6000 7000 

Wavelength (A) 



8000 



9000 



Fig. 1. — Example spectra of a solar-metallicity M d warf (dM-red ) , a su bdwarf (sdM-orange), extreme subdwarf (esdM-green) and 
ultra-subdwarf (usdM-blue). The boundaries between the Lcpinc ct al. (2007) metallicity subclasses are given in the legend. Major atomic 
and molecular spectral features are marked. The CaH near 6800 A grows stronger with decreasing metallicity, while TiO near 7100 A 
weakens. The Na I doublet near 8200 A also weakens with declining metallicity. 




1 1 I I 

sdM + esdM + usdM ■ 
Bochanski et al., 201 1 
Bochanski et al., 2010 ■ 




Fig. 2. — The g — r,r — z color-color locus of SDSS DR9 stars. 
The stars with clean and precise photometry were chosen, and the 
main sequence was selected by selecting stars with g — r < 1.8. The 
(g — r) dM fit to the median g — r, r — z colors (blue circles) given 
in Equation [4] is over-plotted in green. Increasing values of 5( g _ r ) 
are denoted with the arrow. Low-mass subdwarfs generally have 
redder g — r colors at the same r — z color (or spectral type) and 
would be found above th e green line . This locus also traces the 
C ~ 1 bins in Figure 10 of I West et all H2011bl) . to < 0.05 mag. 



Fig. 3. — The Mr, r — z CMR for nearby stars (open cir- 
cles and red line Bochanski ct al. 2010), the statistical parallax 
analysis of SDSS M dwarfs (green filled circles and dashed line 
Bochanski ct al. 2011) and the statistical parallax analysis of SDSS 
subdwarfs (blue filled circles and line). The subdwarf LP 21-194 is 
plotted as the open star. Note the separation in absolute magni- 
tude between the subdwarfs and disk stars at a given r — z color. 
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very little u-band flux. As noted in iWest et all (|2004[) . 
the g — r colors of subdwarfs are redder at a given spec- 
tral type or r — i color, which prompted the construction 
of the S(g-r) index. Both £ and 8r g - r \ are crude metal- 
licity indicators, but they do trace the same behavior, as 
shown in Figure [7] There is a correlation between the 
two values, albeit with large scatter, with metal-poor 
values of C corresponding to larger <5( fl _ r ) values. 

The right panels of Figures [4j [SJ and [6] show the CMRs 
derived as a function of 5r g - r \. The 5( g - r ) CMRs share 
similar trends when compared to the CMRs derived for 
different metallicity subclasses. Stars with larger 5r g ^ r \ 
values are intrinsically fainter, as expected for more 
metal-poor stars. The separation between ^( g _ r ) loci 
is smaller, suggesting a finer sampling of the metallicity- 
absolute magnitude phase space. 

We derived photometric parallax relations by fitting 
the M r i Z ,r—z, <5( ff _ r ) plane, with the following functional 
form: 

M rt i, z = a + ax x (r - z) + a 2 x (r - z) 2 

+a 3 x 5(g- r ) + a,i x (r - z) x 5( ff _ r ), (5) 

where a^.-.a^ are given in Table [4] We selected simple 
functions: low-order polynomials and cross-terms to de- 
rive the relations, while still producing good fits for each 
relation. The standard deviation between the data and 
each fit was ~ 0.1 mag. Combined with the uncertainty 
of each point ~ 0.4 mag, the total uncertainty is 0.41 
mag. The functions, standard deviations and applicable 
color ranges are listed in Table [U These parallax re- 
lations are most appropriate for SDSS observations, but 
should be applicable to any calibrated ugriz photometry. 



4.2. Kinematics of M Subdwarfs 

In addition to estimating the absolute magnitude of a 
homogeneous set of stars, the statistical parallax method 
also measures the velocity ellipsoid (dispersions and ori- 
entation) and the mean motion of the population. The 
mean motions are reported as the Sun's peculiar velocity 
with respect to the statistical parallax targets. We com- 
pared the mean motions of the s ubdwarf sample to disk 
M dwarfs (Bochan ski et al.l 120 1 If ) as a function of r — z 
color. 

4.2.1. Solar Peculiar Motion 

The reflex motion of the Sun is reported by the sta- 
tistical parallax analysis, measured with respect to the 
mean velocity of the subsample. Subsamples with dif- 
ferent mean velocities will reflect different values for the 
Sun's peculiar motion. We employ a right-handed co- 
ordinate system with U increasing towards the Galactic 
center, V increasing in the direction o f solar motion, and 
W in creasing vertically upward (as in iDehnen &: Binnevl 
1998). In this coordinate system the angular momen- 
tum vector of the solar orbital motion points towards 
the South Galactic Pole. 

In the left panel of Figure [8] we compare the mean 
motions of the metallicity subclasses as a function of 
r — z color. The mean solar motion m easured with re- 
spect to disk M dwarfs is also included ([Bochanski et al.l 
120 111) . For most color bins, the motion of the entire 



subdwarf sample (sdM + esdM + usdM) closely mim- 
ics the motions of the sdM population, since most of the 
stars in the sample are found in this subclass. The U 
direction shows a large dispersion among the subclasses 
(i.e., r — z = 1.05), but with general agreement between 
the subdwarf and main-sequence populations, indicating 
similar radial motions through the Galaxy. The W di- 
rection exhibits large dispersions within individual color 
bins as well, but the entire subdwarf population shows 
agreement with the disk M dwarfs at the ~ 5 km s _1 
level. 

The most significant contrast between the disk M 
dwarfs and subdwarfs is manifested in the V direc- 
tion, where the Sun's reflex motion is near 200 km 
s _1 for the bluest subdwarf subsamples. Past inves- 
tigations of subdwarf motions (i.e., iGizisI 17997; have 
measured ( V ) ~ 200 km s _1 . This is very simi- 
lar to the Sun's cir cular velocity around the Galactic 
center (220km s" 1 ; iKerr fc Lvnden-Belll fl986h . Thus, 
the subdwarfs comprise a population with very little or- 
dered rotation, indicative of halo orbits. These orbits 
are characterized by high eccentricity and inclination. 
Orbits for ultra-cool subdwarfs (spectra l types > M7) 
have been a recent topic of investigation (|Cushing et al.l 
[2001 IBurgasser et al]l200lll2009h . since kinematics may 
be useful for constraining ages of brown dwarfs, or iden - 
tifying new tidal debris streams (Bclokur ov et al.1 12006). 
At redder colors (and closer distances) , the average V ve- 
locity of the subdwarf population falls to ~ 60 km s _1 . 
This velocity is significantly larger tha n nearby M dwarfs 
(24 ±3 kms" 1 ; iBochanski et al.ll2011[ ). but smaller than 
the bluest subdwarfs. These stars are probably members 
of the thick disk, while the bluer subdwarfs are primarily 
halo constituents. The mean velocities of each subsample 
can be found in Table [5] 

4.2.2. Velocity Dispersions 

Subdwarfs comprise a kinematically hotter population 
than disk M dwarfs, as shown by their larger velocity 
dispersions. In the right panels of Figure [U we compare 
the subdwarf metallici ty classes to the M dwarf sample 
([Bochanski et al.ll2011[ ) as a function of r — z color. Some 
immediate trends are obvious. First, the total subdwarf 
sample (open diamonds) has larger velocity dispersions 
at all color bins, when compared to the M dwarf sam- 
ple (red filled circles). While velocity dispersions within 
the disk are thou ght to increase roughly as the square- 
root of age (e.g.. iFuchs et al.l 120011 : iHanninen fc Flvnnl 
I2002D due to disk heating, the halo dispersion is gov- 
erned by the Galactic potential. This implies that veloc- 
ity dispersions do not closely follow the relative age dif- 
feren ces. However, both observatio nal dMonteiro et al.1 
2006T) and theoret ical investigations (Zolo tov et al.ll2009t 
Brook et al.ll2012h demonstrate that halo subdwarfs com- 
prise an older, metal-poor population than nearby disk 
stars. The bluer r — z bins contain subdwarfs with dis- 
persions of (err/, cry, aw) ~ (110,100,70) km s~ x . The 
typical uncertainty in each dispersion measurement is 
^5 — 10 km s _1 . These disp ersions suggest the stars are 
drawn from the inner halo (|Carollo et al.|[20Tol) . More- 
over, the radial (U) and vertical (W) dispersions in- 
crease with decreasing metallicity f or most bins, follow- 
ing trends seem in other studies ([Carollo et al.l 120101 : 



TABLE 4 

Absolute Magnitudes of SDSS M Dwarfs 



Absolute Mag. 






Coefficients 51 






^Tot b 




Color Range 






ao 


ai 




as 


a 4 










M r 

M z 


7.9547 ± 1.578 
7.3478 ±1.503 
7.3798 ± 1.6723 


1.8102 ± 2.294 
2.0068 ± 2.189 
1.4977 ± 2.430 


-0.17347±0.814 
-0.45371± 0.776 
-0.33956 ± 0.862 


7.7038 ± 4.511 
7.0445 ± 4.249 
8.4248 ±4.738 


-1.4170 ± 2.771 
-1.2128 ±2.612 
-1.9804 ±2.908 


0.41 
0.41 
0.42 


1.0 < r - 
1.0 < r - 
1.0 < r - 


z < 2.0; 0.0 < o (g _ r) 
z < 2.0; 0.0 < 5( s -r) 
z < 2.0; 0.0 < <5 (9 _ r) 


< 0.5 

< 0.5 

< 0.5 



Note. — The typical uncertainty for the reported absolute magnitudes is o~m — 0.42, which is directly computed from Equation fll 
a The functional form for the photometric parallax is given in Equation [5] 
k The uncertainty is given by CTot — \f a f\t "^^"m' wnerc a M ~ 0.4 and a^ lt ~ 0.1. 
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Fig, 4.— Left Panel: M r vs. r — z for vari ous C cuts. The open circles are nearby stars with accurate trigonometric parallaxes from 
Bochanski (2008) with the best-fit relation of Bochanski ct al. (2010) shown in red. The subdwarf LP 21-194 is plotted as the open star. 
The solid black line denotes the entire subdwarf sample, with various £ cuts described by the legend. The spectral metallicity classes are 
well behaved, with clear separation for most color bins. Right Panel: M r vs. r — z as a function of Sr g _ r \, with same symbols and lines as in 
the spectral type panel. The separation between adjacent 5( g _ r ) CMRs is smaller, indicating sampling of metallicity-absolute magnitude 
space. 




Fig. 5. — Left Panel: M; vs. r — z for various f cuts. The lines and symbols follow Figure[4] Right Panel: M; vs. r — z as a function of 
5/g— r ), with same symbols and lines as in the spectral type panel. 
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Fig. 6. — Left Panel: M z vs. r — z for various f cuts. The lines and symbols follow Figure[4] Right Panel: M z vs. r — z as a function of 
S(g_ r \, with same symbols and lines as in the spectral type panel. 
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Fig 



7. — f vs. (5( g _ r ) for the subdwarf sample. Note that 5( g _ r ) 
roughly traces with small values of £ corresponding to larger 
(redder) values of <5( g _ r ). Each point is color— coded by the star's 
r — z color according to the color bar. 



IBond et al.M2010h . Similar to the M dwarfs, the smallest 
velocity dispersions are seen in the reddest three r — z 
bins in the W direction, indicating that subdwarfs with 
these colors may possess orbits with some ordered rota- 
tion, and suggesting that these stars are drawn from the 
thick disk. In Tabled we list the velocity dispersions for 
each spectral subclass as a function of r — z color, as well 



(2011). 



4.2.3. Probability Plots 



The statistical parallax method assumes that each 
subsample is composed of a kinematically homogenous 
population which can be described by a velocity el- 
lipsoid. Some metal-poor stellar samples have con- 
tained stars from the thick disk and halo, which possess 
different unde rlying kinematics (jBochanski et alJl2007at 
ICarollo et all l2010h . If a given subsample was com- 
prised of similar numbers of thick disk and halo subd- 
warfs, a simple one component fi t would not be valid . 
We constructed probability plots (jLutz fc UpgrenlH98Ct 
iReid et aT1l2002f ) to test the validity of this assumption. 
Space motions were calculated from the proper motions 
and radial velocities from SDSS, along with distances 
estimated from the M r relation in Equation [5j An il- 
lustrative example of the probability plots are shown 
in Figure [9l The U and W probability plots for the 
r — z = 1.05 bin for the various metallicity subclasses 
are plotted. Probability plots are related to cumulative 
distribution functions, with the abscissae calculated in 
units of standard deviation. Hence, a Gaussian distribu- 
tion would be plotted as a line, with a slope equal to the 
dispersion. Populations consisting of multiple Gaussians 
would be manifested as a segmented line (see Figure 3 of 
iBochanski et al.l 12007a! for an example) . The probabil- 
ity plots show that the underlying distributions are well 
described by a single component, validating our assump- 
tion. 

5. CONCLUSIONS 

We have described the statistical parallax analysis 
of low-mass subdwar fs contained in the SDSS catalog 
(|Savcheva et al J 120131) . We divided our sample by color, 
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Fig. 8. — Left Panels: The mean velocities of each metallicity subclass as a function of r — z color. The total sample is marked by 
the open diamonds, and disk M dwarfs are shown with the red filled circles. Each metallicity class is described in the legend. While the 
disk M dwarfs have a smaller mean V velocity (middle panel), the bluest subdwarf bin has V ~ 200 km s _1 , indicative of a non-rotating 
population. Right Panels: Velocity dispersions along the UVW axes for each metallicity subclass as a function of r — z color. Symbols 
are the same as the left panels. Note the larger velocity dispersions at all color bins for the subdwarf population, indicative of an older, 
dynamically heated population. 

metallicity subclass, and photometric properties to fur- 
ther explore the kinematics and intrinsic brightnesses 
of these metal-poor stars. We developed a photomet- 
ric proxy of metallicity for low-mass subdwarfs, derived 
photometric parallax relations, and explored the mean 
motions and dispersions with respect to their metal-rich 
counterparts. 

Metallicity can influence the color-magnitude relations 
of low-mass stars. While this has been directly mea- 
sured fo r more massive stars in the uqr iz photometric 
system (|An et al.ll2008t llvezic et al.ll2008[ ). there has not 
been a detailed study for low-mass stars. We ex panded 
on the initial analysis of iBochanski et all ()201lD to in- 
clude subdwarfs with lower metallicity, as measured by 
the proxy £. We demonstrated that different metallicity 
subclasses (sdM, esdM and usdM) cleanly separated in 
various color-magnitude diagrams. The one nearby sub- 
dwarf with a precise trigonometric parallax and ugriz 
photometry agrees with our derived relations. 

While deriving color-absolute magnitude relations for 
various spectral subclasses does demonstrate the influ- 
ence of metallicity on the intrinsic brightness of stars, 
it requires that spectroscopic observations are available. 
The largest current and future surveys, such as PanStarrs 
and LSST, will be purely photometric, limiting the ap- 
plication of the spectral type relations. We ameliorated 
this situation by deriving 8r g - r ), a photometric proxy 
for metallicity based on a star's g — r and r — z colors. 
Low-mass subdwarfs have red der q — r colors th an their 
solar-metallicity counterparts (jWest et al.ll2004T ) and we 
quantified this effect by measuring the g — r,r — z color- 
color locus for main sequence M dwarfs. The photomet- 



to 
E 




Fig. 9. — Probability plots for the U and W velocity distributions 
of the r — z = 1.05 bin. Each line is described in the legend. The 
various metallicity subclasses are well described by a single compo- 
nent, validating our assumption of kinematic homogeneity within 
each subsample. The V velocity probability plot was omitted since 
the distribution is non-Gaussian due to asymmetric drift, an effect 
which is included in the statistical parallax analysis. 
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ric metallicity parameter <5( g _ r ), reasonably traces the 
spectroscopically determined ( values for low-mass sub- 
dwarfs, with so me scatter, and co mplements earlier in- 
vestigations by iWest et"afl (|2011b[ ) for solar metallicity 
M dwarfs. These relations have potential to be a powerful 
tool in photometrically identifying low-mass subdwarfs. 

The kinematic properties of subdwarfs were also ex- 
plored in this study. The entire sample was compared 
to disk M dwarfs, and some clear differences were seen. 
First, the mean V velocity of bluer subdwarfs is near 
200 km s , suggesting this population is largely non- 
rotating. These objects possess halo-type orbits, with 
high eccentricities and inclinations. The velocity disper- 
sions of the bluest (in r — z) subdwarfs was similar to 
other measurements of the Milky Way's inner halo. The 
redder r — z bins contained stars with kinematics similar 
to the thick disk. 

The Milky Way's halo has yet to be mapped by its 
most numerous inhabitant, low-mass subdwarfs. While 
existing surveys, such as SDSS and 2MASS have offered a 
small glimpse into the low-mass components of the halo, 
this phase space is largely unexplored. The relations de- 
rived in this study will enable easier identification and 
characterization of these objects in the next generation 
of astronomical surveys, such as LSST and GAIA. These 
surveys will probe the halo in situ and offer the most 
detailed picture of its kinematic and structural composi- 
tion. 
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TABLE 5 
Mean Motions of M subdwarfs 



£7 (kn 



V (kn 



W (km s" 1 ) 





sdM 


csdM 


usdM 


Total 


Disk 


sdM 


csdM 


usdM 


Total 


Disk 


sdM 


csdM 


usdM 


Total 


Disk 


1.05 


-14 ± 21 


21 ± 26 


43 ± 1 


19 ± 12 


8 


179 ± 23 


196 ± 23 


187 ± 1 


180 ± 15 


30 


3 ± 11 


15 ± 10 


2 ± 1 


8 ± 8 


5 


1.35 


6 ± 9 


16 ± 9 


20 ± 1 


11 ± 1 


8 


102 ± 12 


180 ± 17 


144 ± 1 


132 ± 1 


25 


8 ± 5 


-1 ± 1 


-18 ± 1 


1 ± 1 


4 


1.65 


6 ± 4 


13 ± 2 




5 ± 3 


8 


53 ± 5 


170 ± 43 




69 ± 5 


23 


3 ± 3 


-1 ± 1 




2 ± 3 


5 


1.95 


4 ± 4 






5 ± 4 


7 


37 ±5 






44 ± 6 


22 


1 ± 3 






1 ± 7 


7 


2.25 


8 ± 10 






± 44 


8 


51± 14 






59 ± 16 


23 


9 ± 9 






12 ± 9 


7 



No te. — We employ a right — handed coordinate system with U increasing 
(as in lDehnen fc B inncv 1993). 



towards the Galactic center, V increasing in the direction of solar motion, and W increasing vertically upward 



O 

B 

a 

c 

CO 

a 

CO 
CO 



■CO' 

f 
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TABLE 6 

Velocity Dispersions of M subdwarfs 



r — z 




ay (km s 


- 1 ) 






°V 


(km s 










(km s 








sdM 


csdM usdM 


Total 


Disk 


sdM 


csdM 


usdM 


Total 


Disk 


sdM 


csdM 


usdM 


Total 


Disk 


1.05 


127 


165 114 


136 


48 


76 


104 


131 


110 


64 


63 


59 


90 


77 


31 


1.35 


100 


140 131 


120 


44 


118 


114 


117 


124 


57 


63 


71 


77 


72 


29 


1.65 


68 


130 


80 


47 


92 


100 




103 


45 


47 


77 




53 


28 


1.95 


61 




68 


48 


82 






95 


44 


45 






47 


29 


2.25 


69 




74 


44 


84 






114 


34 


48 






51 


26 



Note. — The typical uncertainty for the dispersions is 5-10 km s 



